Mechanical restriction versus human overreaction triggering congested traffic states 
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A new cellular automaton (CA) traffic model is presented. The focus is on mechanical restrictions 
of vehicles realized by limited acceleration and deceleration capabilities. These features are incor- 
porated into the model in order to construct the condition of collision-free movement. The strict 
collision-free criterion imposed by the mechanical restrictions is softened in certain traffic situations, 
reflecting human overreaction. It is shown that the present model reliably reproduces most empirical 
findings including synchronized flow, the so-called pinch effect, and the time-headway distribution 
of free flow. The findings suggest that many free flow phenomena can be attributed to the platoon 
formation of vehicles (platoon effect). 

PACS numbers: 89.40.-a, 45.70.Vn, 05.45.-a, 05.20.Dd 



Traffic flow phenomena have been analyzed and mod- 
elled from the viewpoint of statistical physics since the 
early 1990s. The main issues are the characteriza- 
tion of traffic phases and transitions among them as 
well as the development of appropriate traffic models. 
The investigation of empirical data has lead to the 
identification of three traffic phases (free, synchronized, 
jammed) 0, 0, El EL IE IS • Furthermore, time-headway 
distributions extracted from single vehicle data have been 
reported || , showing that even free flow is not as trivial 
as previously believed. In order to reproduce the em- 
pirical findings, many traffic models have been proposed 
both from the microscopic |E 0, El an d the macro- 
scopic viewpoints. Also, there are efforts to estab- 
lish a link between them [13, [lj, lla m order to probe 
any universal features. These achievements have deeply 
influenced the understanding of traffic flow far from equi- 
librium with interacting particles showing collective be- 
haviors. 

A primary objective of traffic models is to explain the 
synchronized flow phase that is characterized by a consid- 
erably high flux without any clear density-flux relation. 
Unfortunately, most models produce merely the free flow 
and jammed traffic. The fluctuating profiles, appearing 
during the transient process heading towards the phase 
separation of the two phases above, are regarded as the 
synchronized flow in those models. However, such a be- 
havior is not compatible with the stability of the syn- 
chronized flow found in empirical data. The approach of 
Knospe et al. 01 produces a different type of solution for 
the synchronized flow, which is not a part of any transient 
process. However, the single vehicle data show unrealis- 
tic velocity fluctuations in this solution. Recently, Kerner 
et al. have claimed to be able to reproduce the synchro- 
nized flow |l8j . However, in their model, a phenomeno- 
logical requirement for synchronized flow is directly im- 
plemented, namely the driver's attempt to adopt the ve- 
locity of the vehicle in front within the so-called synchro- 
nization distance. Another shortcoming of existing traffic 



models is that even the free flow is not reproduced so suc- 
cessfully. In 8] , it is reported that a small time-headway 
below 1 sec is frequently observed especially in the free 
flow. Simple tuning of the spatiotemporal scale to achieve 
such a small time-headway results in an extraordinarily 
high flux. A method to bypass this problem is to utilize 
the gap between the next two vehicles in front , lead- 
ing to a time-headway definitely below 1 sec. However, 
the coupling of two vehicles is inevitable in this approach, 
which is not observed empirically [lflj . In this Letter, a 
new Cellular Automaton (CA) traffic model is presented, 
which overcomes the problems stated above. The me- 
chanical restriction is introduced in the form of limited 
acceleration and deceleration capabilities. Also, human 
behavior is modelled as the driver's excessive response 
to the local traffic conditions. The former is a natural 
consequence of given physical limitations while the latter 
is incorporated into the model to reflect human overre- 
action. In the following, all variables are assumed to be 
integer numbers unless stated otherwise. 

First, we introduce the limited capability of accelera- 
tion (a) and deceleration (D). For simplicity, these are 
both assumed to be constant. The idea of limited accel- 
eration is implemented in most existing CA traffic mod- 
els. However, the deceleration limitation has not been 
enforced strictly or systematically. It is important to 
note that the bounded braking capability changes the 
collision-free mechanism entirely. Most CA traffic mod- 
els impose a collision-free condition explicitly by assign- 
ing arbitrary deceleration values required to prevent col- 
lision. However, it is rather natural to view the collision- 
free flow as a consequence of moderate driving instead of 
infinite braking capabilities. For such a physically realiz- 
able flow, we first design a heuristic collision-free driving 
dynamics strictly observing the limited deceleration. The 
starting point of the new CA model is an inequality which 
guarantees safe driving. Here a vehicle prepares for the 
worst case, namely that the leading vehicle may brake 
suddenly at any time t. Since the follower's reaction is 
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delayed due to the response time which is assumed to 
be the unit time of the model, whether secure driving is 
possible or not is determined at time t + 1 when the re- 
action begins. The velocity which allows safe movement, 
represented by c^ +1 below, should satisfy 

r t (4 +1 ) n(vi +1 ) 

4+A+ £ (c t n +1 -Di)<x t n+1 + J2 «+i- Di ), 

i=0 i=l 

(1) 

where (v^) is the location (velocity) of the n-th ve- 
hicle at time t and the increased index represents the 
vehicle in front. A is the minimal coordinate difference 
required by the follower to guarantee its safety, and thus 
assumed to be at least the length of a vehicle L. Each 
summation accounts for successive decelerations during 
time steps i = 0, 1, .., T{ (i = 1, .., t\) with maximum brak- 
ing capability D, where Tf (t\) for the follower (leader) 
will be specified below. The zero-based summation in- 
dex stands for the response time of the follower. For 
T f.i( v ) — V /D an d A = L, Eq. Q suggests such c^ +1 
that can guarantee a complete stop showing bumper-to- 
bumper configuration. We call this dynamics the strict 
collision-free dynamics (or criterion) in this work. Later 
on, the expressions for T{a and A will be modified to 
take into account the human overreaction. Note that for 
a given n, the safe velocity c^ +1 is not unique but just 
has an upper bound. In the following, the largest c^ +1 
satisfying Eq. (Q, denoted by c^ +1 , is used to reflect the 
desire of drivers to move as fast as possible. 

Next, an element of human behavior is introduced 
which is actually responsive to the local traffic situation. 
Generally, it is accepted that a driver's maneuvers are 
not precisely predictable by a simple rule. The usual 
solution for this problem in simulation models is to in- 
troduce fluctuations which cover the human factor in a 
stochastic way. However, a different strategy is adopted 
here. It is supposed that the driver's behavior may be 
biased depending on the local traffic situation. To imple- 
ment this simply, a 2-state variable is introduced: 

t = | for < v* n+1 < v l n+2 or v* n+2 > v fast , ,^ 
1 1 otherwise 

with a constant Uf ast slightly below v ma x- The state 7^j = 
corresponds to a situation where the driver judges that 
the local situation is optimistic since the cars in front are 
speeding away. It is assumed that the driver will move 
faster in order to catch the car ahead, even faster than the 
allowed velocity under the strict collision- free dynamics. 
This state is denoted as the optimistic state. Otherwise, 
for 7„ = 1, the driver is in the defensive state, where the 
vehicle slows down below the velocity due to the strict 
collision-free criterion. This distinction in respect to the 
local traffic situation is termed human overreaction here. 
It is realized in the model by manipulating T{ t i(v) and A 



in as follows: 

A = L + 7 r t l max{0, min{g add , v l n - # a dd}}, 
T f {v) = J* n v/D+ (1 -7*)max{0,min{u/AWe} - 
n (v) = 7* v/D + (I - 7* )mm{u/.D, i safo }. 

(3) 

Herein g ad d is introduced for an additional security gap 
in the defensive state (7* = 1) and t sa f c is a maximal 
time step during which the follower observes his/her own 
safety in the optimistic state. The additional —1 for T{(v) 
compensates for the surplus time step due to the fol- 
lower's response time only when 7^ = 0, and thus the 
role of i sa f e is properly implemented. For 7^ = 0, the Tfj 
can be smaller than those (v/D) necessary for complete 
stops while A returns to L. Consequently, lower safety is 
required compared to that needed by strict collision-free 
dynamics, and thus a faster c^ +1 can be chosen. On the 
other hand, for 7^ = 1, A can be larger than L while 
Tf j i(u) return to v/D, which implies over safety. In this 
way, a lower c^ +1 can be assigned. 

The update rules of the model can be written in the 
following form: 

1? = max{p d ,po - Vn(P0 -Pd)/«slow} 

2 &+ l = max{4 +1 | 4 +1 satisfies Eqs. (HH3J} 

3 = min{w max , v l n + a, max{0, v l n - D, c^ 1 }} 

4 = max{0, v l n - D, v 1 ^ 1 - i]} 

where rj = 1 if rand() < p, or otherwise 

C ™t + l n*t _l_ 

Herein, the stochastic parameter p (< I) in step I lin- 
early interpolates between pg and pd if v ^ is smaller than 
f4iow (note that p is a real number). We set po > pd so 
that step I is a g eneralization of the well known slow-to- 
start rule [171 |2C1 121| , which is known to be an ingredient 
in the formation of congested traffic states. Step 3 guar- 
antees that the updating velocity satisfies the mechanical 
restriction as well as the traffic regulation. The stochas- 
tic deceleration in step 4 is also limited by the braking 
capability D. Note that steps 3 and 4 observe the lim- 
ited deceleration consistently • The length of one cell 
is chosen to be Aa; = 1.5 m and the unit time is set to 
At = 1 sec. The following model parameters are moti- 
vated by empirical facts: a — 1, D = 2, L = 5, Wf ast = 19, 
tsafc = 3, g ad d = 4, po = 0.32, p d = 0.11, Usiow = 5, and 

%ax — 20. 

In Fig. ^ the traffic phases occurring under periodic 
boundary conditions are depicted. Three traffic states 
(free flow, synchronized traffic, jams) can be identified in 
the fundamental diagram Fig. a) (see 2]). The straight 
line with the positive slope corresponds to free flow. The 
synchronized states form a 2-dimensional region in the 
middle of the diagram while jammed vehicles produce the 
scattered points below. A typical spatiotemporal shape 
of each phase is shown in Fig. GJb) . Note that the so- 
called universal constants of traffic flow 0, 0, 0] , such as 
the jam velocity v g w —15 km/h and the flux out of a 
jam q out ~ 1800 veh/h, are also reproduced by simply 
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FIG. 1: Results for a road (periodic boundary conditions) con- 
sisting of 40,000 cells (60 km). A homogeneous distribution 
of standing vehicles is used as the initial condition while the 
density is varied from 10-50 veh/km with a step of 2 veh/km. 
(a) Fundamental diagram: Single vehicle data is gathered at a 
fixed position and then the flux (J) and velocity (v) are aver- 
aged every minute. The density (p) is simply obtained via the 
hydrodynamic relation J = pv. The measuring time is 10, 000 
sec after the relaxation of 30, 000 sec. (b) The spatiotemporal 
shape of the different traffic phases [F(ree), S(ynchronized), 
J(ammed)] is depicted for 30 min. The initial densities are 
16, 30, and 44 veh/km, respectively, (c) Snapshot of synchro- 
nized flow (randomly chosen part zoomed at 40, 000 sec) with 
an initial density of 30 veh/km. The filled squares represent 
vehicles moving from left to right. 

adjusting p and pd- For the purpose of demonstrating 
the inner structure of synchronized flow a snapshot of a 
part of the road is presented in Fig. ^c). As required, 
this is not a transient process and exhibits a smooth ve- 
locity profile. This implies that the points comprising the 
synchronized area in Fig. ^a) are not attributed to the 
averaged effects of strong fluctuations but to the special 
headway-velocity relation. 

The results for an open system with on-ramp are pre- 
sented in Fig. [3 We first examine the so-called pinch 
effect 0, 0, Il8| describing the process of a local self- 
compression in synchronized regions, which leads to the 
formation of small narrow jams. These small jams evolve 
finally into a few wide jams through a merging process 
while moving upstream when the small jams grow enough 
to lose the stability of the synchronized flow. This pro- 
cess can be seen in Fig.EJa). Empirically, such a process 
is the most frequent type of congested traffic near a bot- 
tleneck and is thus named the general pattern (GP) 0. 
This process is captured more clearly in Fig.^b). Start- 
ing from the top, it is seen that synchronized flow (high 
flux and slow velocity) is formed near the on-ramp. In 
the following two pictures, many small drops are merged 
into small narrow jams along the upstream. Finally, a few 
wide jams remain far away upstream from the on-ramp. 
The other known types of congested traffic near bot- 
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FIG. 2: The impact of an on-ramp at 20.4 km (black arrow) 
with a 0.6-km long merging road is demonstrated. The left 
(right) end of the road is positioned at —45 km (45 km). The 
influx at the left is kept constant during simulation. After q ln 
is supplied for 10 hours, the simulation time is set to and 
then 8 minutes later the on-ramp injection q on is turned on. 
Within the merging road the widest leading gap is selected 
and then a vehicle is inserted at the mid point. For its ve- 
locity, 0.7 times that of the follower's is assigned. Below, a 
coupled number means (qin,Qon). (a) General pattern (GP) 
at (1800,550). (b) Pinch effect in (a): Each figure shows one 
minute averaged velocity (solid line) and flux (dotted line) 
at different locations. From the top, detectors are located at 
20.4, 18.9, 15.9, and 5.4 km, respectively, (c) Localized syn- 
chronized pattern (LSP) flow at (1300,650). (d) Widening 
synchronized pattern (WSP) flow at (1950,350). (e) Moving 
synchronized pattern (MSP) flow at (2050, 150). (f) Phase 
diagram of the congested traffic patterns. The filled triangles 
correspond to mixed patterns. 

tlenecks are localized synchronized flow patterns (LSP), 
widening synchronized flow patterns (WSP), and moving 
synchronized flow patterns (MSP). It is stressed here that 
these patterns are also reproduced in the present model 
as shown in Figs|2Ic,d,e). Furthermore, we remark that 
the phase diagram for the congest ed p atterns [Fig. EJf )] 
is quite comparable to those in 0] • 

Finally, we examine the time-headway distribution of 
the model. The correspondence between the numerical 
results and its empirical counterpart is quite satisfactory 
as shown in Fig.[3fa). In contrast to the result in 01 > the 
peaks below 1 sec are realized merely with identical vehi- 
cles. A typical spatial configuration of the free flow shown 
in Fig. |3{c) indicates that such peaks are attributed to 
the platoon formation of vehicles (platoon effect). Re- 
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(a) (b) 




(c) platoons in free flow 

|* 1 1* * **| *** f*| ***** I I* |**;**|* m I i *l v I v iTtttt^ I'lrl 1 1 1 If! 1 1*! I & 

O 3 
location (km) 



ular platoons as a non-trivial element of free flow. Thus 
the platoon effect is proposed to be the origin of many 
features of free flow. We remark that some situations re- 
quiring greater deceleration beyond the braking capacity 
for safety, such as careless insertion near on-ramp, can 
result in collisions in this model. In summary it can be 
stressed that the presented model sheds light on open 
questions of traffic modelling and is therefore useful for 
the understanding of certain traffic phenomena as well as 
for applications such as fully automated driving. 

H. K. L. was supported by the Brain Korea 21 Pro- 
gram. The authors are grateful to the BMBF within the 
project DAISY for financial support. 



FIG. 3: (a) [(b)] Time-headway distribution of free [synchro- 
nized] flow is depicted. The cross, filled diamond, and open 
diamond represent densities of 12, 16, and 20 veh/km [24, 
32, and 40 veh/km]. The insets are taken from || (empirical 
data). Solid, dashed, and dotted lines correspond to density 
ranges from ~ 12, 12 ~ 24, and 24 ~ 36 veh/km [24 ~ 36, 
36 ~ 48, and 48 ~ 60 veh/km]. (c) Randomly chosen part of 
the road whose overall data give the connected open diamonds 
in (a). 



visiting Eqs. QJ, J2Jl, and © with the requirement of 
free flow (c^ +1 = — u max for all n) helps us to un- 
derstand the formation of such platoons. Then Eq. Q 
allows even 0.3 sec as a time-headway, which opens up the 
possibility of the formation of vehicular platoons in free 
flow. It is emphasized here that these platoons can ex- 
plain the small time-headway frequently observed in free 
flow without an unrealistic high flux. Thus the platoons 
may be one of the fundamental objects which can charac- 
terize free flow. Another important role of the platoons 
is suggested as follows. Since the large gaps between the 
platoons can absorb fluctuations propagating backwards, 
they stabilize free flow and thus influence the stability 
and transition properties. For the case of synchronized 
flow, qualitative agreement with empirical findings |8( is 
observed as follows [see Fig. EJb)]: i) distribution shifts 
to the left as density increases, ii) peaks are near 1.5 sec, 
iii) time-headway smaller than 1 sec still exists, and iv) 
large time-headway events are reduced compared to those 
of free flows. 

In conclusion, a new CA traffic model focusing on the 
mechanical restriction realized by limited acceleration 
and braking capabilities is introduced. A further element, 
namely human overreaction, is implemented in order to 
reflect the driver's tendency toward biased reaction ac- 
cording to the local traffic conditions. It is shown that 
the model reproduces most empirical findings including 
the three known traffic phases, the so-called pinch effect 
and several types of congested traffic patterns as well 
as the small time-headway below 1 sec especially in free 
flow. In particular, the presented model produces vehic- 
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